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Abstract 

We study the implications of the experimental results on the /i — > decay rate and the 
muon anomalous magnetic moment, on muonic lepton flavor violating processes, such as 
/i — > 3e and fiN — > eN. We use a model independent approach in this analysis, where these 
processes are considered to be loop induced by exchanging spin-1/2 and spin-0 particles. 
We explore two complementary cases, which has no or has an internal (built-in) cancellation 
mechanism in amplitudes. Our main results are as following, (a) Bounds from rates are 
used to constrain parameters, such as coupling constants and masses. These constraints can 
be easily updated by simple scalings, if the experimental situations change, (b) The muon 
g — 2 data favors non-chiral interactions, (c) In fi — > 3e and fiN — > eN processes, Z-penguin 
diagrams may play some role, while box diagrams contributions to fi — > 3e are usually 
highly constrained, (d) In the first case (without any built-in cancellation mechanism), 
using the recent fi — > bound, we find that fx — > 3e and fiN — > eiV rates are usually 
bounded below the present experimental limits by two to three orders of magnitudes in 
general. Furthermore, by comparing Aa M and £>(/i — > ey) data, the couplings of ji and e 
are found to be highly hierarchical. Additional suppression mechanism should be called for. 
(e) In the second case (with a built-in cancellation mechanism), mixing angles can provide 
additional suppression factors to satisfy the Aa M and B(fi — > e^y) bounds. While the fi — > 3e 
rate remains suppressed, the bounds on fiN — > eN rates, implied from the latest // — > 
bound, can be relaxed significantly and can be just below the present experimental limits. 



I. INTRODUCTION 



Charge lepton flavor violating (LFV) processes are prohibited in the Standard Model (SM) and, 
hence, are excellent probes of New Physics (NP). Recently the search of /i — > cy decay was reported 
by MEG collaboration giving pQ 

B(^ + -»• e+ 7 ) < 2.4 x 1(T 12 . (1) 

The bound is several times lower than the previous one PJ. This result received a lot attentions (see, 
for example, [3H5]). In many New Physics models this decay mode is closely related to other lepton 
flavor violating processes, such as /i + — > e + e + e~ decays and fj,~ N — > e~N conversions [6]. The 
present limits and future experimental sensitivities [U [21 [7] of these LFV processes are summarized 
in Table [TJ Note that present bounds on [i LFV rates are roughly of similar orders. It will be 
interesting to see what are the implications of the new B(/j, —> ej) bound on these LFV processes 
and the interplay between them. 

Since 2001, the muon anomalous magnetic moment remains as a hint of a NP contribution (see, 
for a review, |5]). Experimental data deviates from the Standard Model (SM) expectation by more 
than 3a |2|: 

Aa /t = a^ xp - a^ M = (287 ± 63 ± 49) ± HT 11 . (2) 

Since NP contributes to Aa^ and £>(/x + — > e + j) through very similar loop diagrams [see Fig. 1(a) 
and (b)], it is useful to compare them at the same time. 

The Large Hadron Collider (LHC) is working well. So far no NP signal is found (see, for 
example [9]). Plenty of well studied NP models or scenarios are ruled out or cornered. Therefore, 
it will be useful to study low energy effect, when the NP scale is still beyond our reach. Given the 
present status on NP models, we believe that it is worthy to use a model independent approach. 

In this work we consider a class of models that muon g — 2 and various muon lepton flavor 
violating processes, such as \x — > ej, fi — > 3e and [i — > e conversions, are loop-induced by exchanging 
spin-1/2 and spin-0 particles. We try to see where the present g — 2 and fi — )■ &y experimental 
results lead us to on estimating rates or bounds on various LFV muonic decay modes and the 
interplay between them. Two cases, which are complementary to each other, are considered. In 
the first case, there is no any built-in cancellation mechanism among amplitudes. The second 



TABLE I: Current experimental upper limits and future sensitivities on various muonic LFV pro- 
cesses piEicrj. 

current limit future sensitivity 

B(fi + -»• e+ 7 ) < 2.4 x 10~ 12 10~ 13 

B(n+ -»• e+e+e") < 1.0 x 10" 12 10~ 14 - 10~ 16 

B(p-Ti -¥ e"Ti) < 4.3 x lO" 12 10~ 18 

B(/i~Au -»• e - Au) < 7 x lO" 13 lO" 14 - lO" 16 

B(pTAl -»■ e~Al) • • • 10~ 16 
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FIG. 1: (a) and (b): Penguin diagrams that contribute to muon g — 2, fi + — > e + 7, /x" 1 
and N — > e~N processes. Note that diagrams involving self energy parts are not shown, (c) 
and (d): Box diagrams contributing to the fi + — > 3e process. Figure (d) takes place only when 
ip m n are Majorana fermions. 



case is with some built-in mechanism, such as Glashow-Iliopoulos-Maiani (GIM) or super GIM 
mechanism. These two cases will be compared. 

The lay out of this work is as following. In the next section, the framework is given. Numerical 
results are presented in Sec. Ill, where bounds from rates are used to constrain parameters, 
such as coupling constants and masses. Correlations between different processes are investigated. 
Discussion and conclusion are given in Sec. IV and V, respectively. Some formulas and additional 
informations are collected in Appendices. 



II. FRAMEWORK 

In this section, we begin with introducing the lagrangian of a generic interaction involving lep- 
tons, exotic spin-1/2 fermions and spin-0 bosons. Formulas of processes of interested, discussions 
on subtleties on the calculation of the Z-penguin amplitude and explicit expressions of Wilson coef- 
ficients will be given subsequently. This section end after the formulation of the two complementary 
cases as briefly mentioned in Sec. I. 
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A. The interacting lagrangian and diagrams 

The lagrangian of a generic interaction involving leptons (I), exotic spin-1/2 fermions (ip n ) and 
spin-0 bosons (fa) is given by 

Ant = MgfiP L + 9&Pr)i$ + MlPr + qirPl)^, (3) 

where summation over indices are understood unless specified. The lagrangian is given in the mass 
bases and are ready to be used in calculations. However, it is important to make sure that it 
transforms as a singlet under the SM gauge transformation. 

In the weak bases of ipLp, ipRp, 4>La and 4>R a , the interacting lagrangian is 

Ant = (g' P iti>R P hfa La + g'm^LplR^Ra) + h.C, (4) 

where we denote <j>L(R) f° r the scalar fields that couple to Iur) and subscripts p and a are the labels 
of different weak fields, which may have different SU(2)xU(l) quantum numbers. It is important 
to require that Ant transforms as a singlet under the SM gauge groups and quantum numbers of 
tp and <fi are related (see Appendix |A|) . 

The mass bases are related to the weak bases through the following transformations 

fa = Uta^La + U%<P Ra , VnL(i?) = V$ R) 1p L (R)p, (5) 

where % and n are labels of mass eigenstates and U and V are the mixing matrices relating weak 
and mass eigenstates. With 

n ni _ /pa v R(L)ttL(R) ( „s 
SlL{R) - 9 i L ( R ) Vnp u ia \°) 

the interacting lagrangian is now brought into the form shown in Eq. ([3]), which is more convenient 
and will be used in later calculations. 

These interactions will induce lepton flavor violating processes at one-loop level. Penguin di- 
agrams contributing to the muon anomalous magnetic moment and /i + — > e + 7 are shown in 
Fig. [l|a) and (b), while box diagrams contributing to the fi + —> 3e process are shown in Fig. [l^c) 
and (d). Note that (i) penguin diagrams shown in Fig. [l|a) and (b) also contribute to fi + — > 3e 
and fiN — > eN process by connecting the virtual photons or the Z bosons to electron currents and 
quark currents, respectively, 1 (ii) Fig. [ljd) takes place only when V'm.n are Majorana fermions. 

B. Formulas for various processes 

To define our variables and to specify our convention, we collect formulas for various precesses 
here. The relevant effective lagrangian in this study is 

£eff = A'Z 7 + Cl'lll + £l>lqq (7) 



It is possible to have box diagrams with the electron line in Fig. [TJc), (d) replaced by a quark one and 
contribute to fiN —> eN conversions in some cases. We will consider this contribution in the discussion 
section. 
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with ?W = e,/x, r and q denoting quarks. Each term will be specified in below. For I' ^ I, we have 

C Vh = 1' l cj^I r F^A l , r + I'ro^IlF^Ar.l + h.c, (8) 

and 

A L r> = A* R , L , A RV = A* L , R . (9) 

Note that in the case of I' = I, we do not need the additional hermitian conjugated terms in Eq. Q. 
These ^4s are from the so-called F2 photonic penguin and their explicit forms will be given later. 
The effective lagrangians for /' — > 31 decays and /' — > I conversions are [6] 

Cl'lil = 9r'lrl(1' 'rIl)(IrIl) + 9l'rlr(1' 'lIr)(JlIr) 

+9r'rrr{1' r1^1r){Ir1^r) + 9l'lll{1' l ^Il){Tl1,j,Il) 
+9r'rll{1'ri' 1 Ir){Il1^l) + 9l'lrr{V lI^IlWrI^Ir) + h.c, (10) 

C nqq = J2[gLv{q)l , Ll' x h + 9Rv{q)l'Rl^lR]qi^q + h.c, (11) 



q=u,d 



where 



2 P Z Z B 

9M'MNO = e- Qi9m'm°no + 9m'm9i n 5no + 9m'mnoi 

9Mv{q) = e.Q\9M'M + ^9m'm(9% l + 9q R ), (12) 

for M, N, 0=L, R with g^, M from the so-called F\ photonic penguin, gfj, M from the Z-penguin 
and 9m'mnn f rom box diagrams. More details and the explicit forms of these Wilson coefficients 
will be given later. Note that although the above £i'i q q is not the most generic one, it contains all 
the relevant parts that are closely related to Cym and Cin T 

We now collect the formulas for various processes originated from the above Lagrangians. Com- 
paring the effective lagrangians of the lepton g — 2 and the electric dipole moment (EDM), 2 

C g -2 = " J^ A «* Iv^lF^, C EDM = ~\d x la^lF^, (13) 

to the generic expressions in Eq. ([8]), the anomalous magnetic moment and EDM of lepton I can 
be readily read off as 

A ai = -^Re(A RL ), d, = 2Im(A RL ), (14) 



respectively. The /' —> Ij and /' — > III decay rates are given by 

r(P - h) = (m t" m f )3 (l-W + 1 W) (is) 



and [6j 



mf, 



\9R'LRL\ |2 , 00 eA R'L 
o r \9R'RLL\ + 31 



m v 



mf, 11, 



1 i" b V ±± ^ 

lo sa ~ — ) 



mf I 



+16Re ( eA RIL gl, LLL \ + 8Re / eA R , L gl, LRR \ 
+L o R, (16) 



2 We use the convention where = + ieQA^ with e 
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respectively. While the I'N — > IN conversion rate ratio is governed by 



Bi'N^eN = CODV , (17) 

^capt 



with 



0J C 



A% r D 



+ 2[2g* LV (u) + g* LV (d)]vW + 2 [g* LV (u) + 2g* LV {d)]V^ 



+ L^R. (18) 

The numerical values of D, V and Co> cap t are taken from |13|. I14j and are collected in Appendix p] 



2 



C. Z-penguin amplitudes 

The calculation of the Z-penguin amplitude is quite complicate and subtle. Some explanations 
are needed. 

The interaction involving a Z boson is given by 

Cf Qt = -lZ{gf L P L + gf R P R )l-^ Lp ^g^ Lp -i, Rq Zg^ Rq 

-igf La (<f>* La d^ La - d^la^La) ~ i 9 l Ra (tpRad^Ra ~ 8^* R J Ra ) + ... (19) 

with 

9X = ■ n 6 n ( T 3 " Sin2 °WQ)X , (20) 

sin t>w cos vw 

for X = Il(r) j 4>L(R)p an d <t>L{R)a m the weak eigenstates. Since C m t transforms as a singlet under 
the SM gauge group, the gj| of various fields are related through 

9i U m ~ 9 %r(l) p-94,a = 0, (21) 



if the corresponding coupling g'f^R) ^ n 

Eq. Q is non- vanishing. 
Although the couplings 9^ L(R) are diagonal in the weak bases, it may have off-diagonal 
terms in the mass bases. In the mass bases, the interacting lagrangian involving a Z boson is given 
by 

+ (22) 

with 

a Z = tjL Z tj]L ,ttR Z tj\R Z =V L(R) Z y\L(R) (2 o) 

y<pij '- J iay<p La '- J aj ' '-'lailfina^aj ' ^i' L ^ R )mn mp iftpL(R)p P n ' \ > 

The one-loop amplitude for I — > I'Z consists of two diagrams shown in Fig. 1 and two additional 
diagrams involving self-energy diagrams with Z attached to external lines. The resulting amplitude 
is given by (neglecting mi, and q 2 ) 

iM = i^'^* p *+^^ 

„Z X _ _ X „Z X „Z x \D_1/„»)D. i „ n JD_V 



+(9i R SijS mn - g^ L m Jij - gjtijSm^PRligflPL + 5/k-Pr> 

1 / 2 47T \ . o n 2 2 71 ,r2 



+ Y 67T 2 ( - 9 t R mn ~ 9$ Lmn )u'(9VtPR+9Vk*PL) f^L ~ 9?rPr) u Gz {™% m , m % n , ) , 

(24) 
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where d is the number of the space-time dimension, is the Eular number and Fz and Gz are 
the loop functions whose explicit forms shown in Appendix [Bj Note that M is an arbitrary mass 
parameter introduced to balance dimension. We will return to it later. 

Note that the divergent part contained in the first term is indeed vanishing by requiring the 
interaction C- mt in Eq. Q be invariant under the SM gauge group, i.e. we have 

mi* f Z r s: Z x Z sr \ nj /pa* r Z Z Z \ /pa p, 

9VL{R)\9l L(R) °ij°mn ~ 9^ R(L) mn°ij ~ 9(j>ij0mn)9 lL ( R) - 9 v L(R)\9l L(R) ~ 9if} R{L) p ~ 9<f>a)9 l L ( R ) ~ U, 



R(L)P 



(25) 



where sum over indices are understood and Eq. (21) has been used in the last step. The non- 



divergent part, namely the one with Fz, survives. Note that by the same token the dependence 
on the arbitrary mass parameter M cancels. The resulting Wilson coefficients will be given later. 
It is easy to see that in the non- mixing case {U = V = 1), the whole first term (including the Fz 
term) is vanishing. 

As a cross check, we note that the same expression of iM can be use to obtain the lowest order 7- 
penguin amplitude by replacing each g z by the corresponding eQ. Since under these replacement 



tj^mn 



and (g, 



e(Qip — Qip)S mn = 0, the corresponding u' fu term is vanishing as expected. 



i/j^mn 



9tj) R mn 



D. Wilson coefficients 



Induced by the interaction given in Eq. ([3]) the Wilson coefficients for the effective lagrangian 
in the Sec. II A are calculated to be 

Am ' N = ^2^ mi ^M9?M + rnvg^g^iQ^F^ml^m 2 ^) - F x (m^ , ) ) 

+%M(Q^3(mi,m^) - Q^F 2 {m\^m\ n ))\, (26) 

for M 7^ N, but with M, N = L,R, and Fi are loop functions collected in Appendix [Bj The Wilson 



coefficients in Eq. (12) are 



p 
9r'r 



9wr 



9vl 



+m^ n (rmg^gfl + m v g^9m) [Q^ n G 3 (m% , m\ n ) + G 3 (rn\ n , m\, )] } 
~ 167r 2 m| sin 28 w 2KRi i mn9 ™R '9?R F z( m l m > ™\ n , m % > ™% , m z) 
2AT 3i>mn g^*g^G z (m% m , m\ n , m\. ) , 



16ir 2 m z s'm28\Y 
1 



o mi* njrn/ 2 2 2 2 2 \ 

16ir 2 m%sm2d w Lljmn9l ' L 9 iL F z( m ^ m ^ m 4>o m 4>p m z) 



"z • 
e 



16-7r 2 m^ s'm29\Y 



2AT 3ll!mn g™*9?£Gz (m% m , m% n ,m%), 



jj 

9r'rrr 



9r'lrl - Yfr^ ^ m ^m' m Vn' m 0i' m 0i^^'« 9il 9ir 9il ~ ir )9vR 9m 9il9\l 



1 



16vr 2 
1 

~4 



9?r9ir gmgmFim^ , , m . , m^. ) 



mi* mi nj* ni r~H 2 2 2 2 

9vr 9ir9ir giRG(m^ m , , m^. , m 
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with 



9r'rll = ^{ " ^(m^m^X^my^ 

+\9?Rgn*9lR9iL G ( m % m , < . < - <• ) } , (27) 

K L(R)ijmn = Sm29w {9l L{R) ^ij^mn ~ 9^} R{L) mn^ij ~ d^ij^mn) / 'e, 

ATs^mn = V mp T^ R pVp n — V m <pT^ LV Vp n , (28) 

loop functions F and shown in the Appendix |B| and t] = 1(0) for Majorana (Dirac) fermionic 
ip. Other g can be obtained from the above ones by exchanging R and L. Note that for definiteness 
we take M = mz in Fz- As before the summation on m,n,i,j is understood. 

E. Two Cases 

We consider two complementary cases. 
1. Case I 

In the first case, namely case I, there is no built-in cancellation mechanism. The amplitudes 
may contain N different sub-amplitudes, each comes from one of the loop diagrams as shown in 
Fig.0 

N 

A = Y,Aj. (29) 

3=1 

We will constrain parameters from data by switching various diagrams (sub-amplitudes) on one at 
a time. The corresponding Wilson coefficients of a typical sub-amplitude can be obtained by using 



formulas in Sec. II D but with the replacement, 

9i1a -»• 91M, (30) 
with all summation on n and I suspended. Since there is no built-in cancellation in this case, 
different sub-amplitudes are in principle independent from each other. Although it is likely to have 
various amplitudes to appear at the same time in a realistic model calculation and to interfere 
with each others, the interference effects only become important if the amplitudes are of similar 
size. Hence, our analysis not only is valid when the sizes are different (hence, constraining the 
most dominant amplitude), but can also provide information on regions where interference may be 
important. 

2. Case II 

In case II, there is a built-in cancellation such as a GIM or a super- GIM mechanism in the NP 



sector. This case is complementary to the previous one. Some of the sub-amplitudes in Eq. (29) 
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are intimately related. We have to group them to allow the cancellation mechanism to do its job 
first. The grouped amplitudes should be viewed as new sub-amplitudes and we will turn them on 
one at a time to constrain their sizes from data. To be specify, we consider 

9im -> 9im = 9imT z m, (31) 

where we have M = L,R and giu is real as the phase is absorbed into T. Note that the matrix T is 
similar and related to the mixing matrix U, but is not identical to it. These T satisfy the following 
relations: 

M i N \ m {f>)MN> L M L N= d d MN- (32) 



A typical expression of Wilson coefficients given in Sec. II D is transformed in the following way: 

d 

9™Mf{™% mli)9 l eN -> rnl-^fim^mfyg^MgeNSjfe 1 *, (33) 
i om 4> 

where is the average mass squared of fa and the mixing angle 5jf e N is defined in the usual way 

to be ma 



-. ( 2MN\ 



The Wilson coefficients in this case can be obtained readily by applying the above replacements to 



the generic formulas collected in Sec. II D 



Note that in the Z penguin amplitude the zeroth and first order terms in the k-l(R)Fz part are 
vanishing. The leading order contribution is at the level of Slr&rl-, which is beyond the accuracy 
of the present analysis and are neglected. 



III. RESULTS 

Numerical results in cases I and II are given in this section. Unless specified explicitly, experi- 
mental inputs are taken from Table 1 and Ref. [2]. 



A. Case I 

In Fig. [2] we show the allowed parameter space for Qcj),ip\g^L{R)\ 2 l m % an d Q^,^Re((7* / j(? M L)/m^ 
constrained by the measured Aa^ with exclusion of \g^L(R)\ 2 ^ \9hL9iir\ > 47r and < 100 GeV. 3 
The latter requirements are to ensure perturbativity and to satisfy the experimental bounds on 
the masses of exotic particles [2j . Bands denoted with <fi or ip are allowed regions obtained through 



3 It is easy to see that |g M L(i?)| 2 > 4-7T and < 100 GeV implies {g^iR) 1 2 l m % > 47r/(100) 2 

GeW 2 and \g ^r^ I > 4-7r/100 GeV -1 , while \g^L{R)\ 2 > 4-7T and < 100 GeV implies 
\g^L(R) I 2 /™ 2 = (l5 M L(fl)| 2 /™ 2 ) x (m^/m^) 2 > 4tt/(100) 2 x (m /m^) 2 GeV" 2 and \g^ L g^R\ 2 /nty > 
(47r/100) x (m^/m^,) GeV -1 . These excluded regions are shown by shaded areas with horizontal or 
inclined boundaries. 
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\Qw Im(^ R *^ L )l / m^CGeV" 1 ) 




FIG. 2: (a) Allowed parameter space for ±Q<^|g M L(_R)| 2 / m2 constrained by Aa M (bands with 
solid or dashed lines) with exclusion of \g l _ l L(R)\ 2 > 4vr and m^^ < 100 GeV (shaded regions with 
dot-dashed lines), (b) Allowed parameter space for ±Q ( j,^Re(gJ lR g fl L) /m^ constrained by Aa^ 
with exclusion of {g^g^M > 47r and m^^ > 100 GeV. Excluded parameter space (shaded regions 
with solid or dashed lines) of \Q^lxn(g* R g fl ] J )\/m^ from the muon EDM bound and the expected 
sensitivity are also shown. 

contributions from diagrams with <j> or interacting with a photon [see Figga) and (b)]. Excluded 
parameter space of \Q ( f >t1 pIm(gJ lR g fi L)\/m 1 p confronting the muon electric dipole moment (EDM) 
bound [2] is also shown in Fig. [2^b). 

From Fig. j2^a), we see that the allowed regions on Q <f>\g nL(R)\ 2 an d ~Q^\g^L{R)\ 2 / m2 b are 
similar and the signs of Q<f>^ are constrained by data. Note that the allowed parameter space is 
quite limited. Indeed, it is almost closed by the bounds from \g^L(R)\ 2 < 47r and m$ > 100 GeV 
and > 100 GeV. The allowed region is around iQ^lg^i?)! 2 / 771 ^ ~ 10~ 4 ~ 10~ 3 GeV -2 and 
m^/m^ — 0.3 ~ 3, which implies that for of a few hundred GeV, the couplings Q^I^l^)! 2 
are required to be of order 0(1) ~ O(10), which are rather large, and are even larger for heavier 
m^. To see it in another way, if we take the size of g^L^R) to be similar to that of the electric 
coupling e, we need to be as light as 10 to 30 GeV to reproduce the experimental result on 
Aa^. Thus, it is unlikely to use a chiral type interaction (g^L x g^R = 0) to generate the measured 

From Fig. [2^b), we see that the allowed parameter space is substantially larger. To reproduce 
the measured Aa^, the mass ratio has to be in the range of O(10~ 5 ) < m^/m^ < C(10 2,3 ), which 
is much wider than the one in Fig. [2^a). We note that the bands of the allowed parameter space 
behave rather differently in two regions roughly separated by m^/m^ = 0.1. (i) For m^/m^ ^$ 0.1, 
the horizontal bands denoting the allowed parameter region for ±Q t p^Ke(g* R g IJi i J )/m^ are around 
4 x 10 -6 GeV -1 . They are insensitive to m^/m^, since the chiral enhancement factor m^/m^ 
compensates the suppression from the heavy ip mass. Note that m^/m^p can be as low as 3 x 10~ 5 , 
which implies that up to 3 x 10 3 TeV is still capable to reproduce the measured Aa M in the 
extreme case, where we have a light <f) [m^ = 0(100) GeV] and large couplings [\g^Lg^R\ = 0(4tt)]. 
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. ""7i->3e, 




^->ey, 
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(d) 



f~i i<r j - 

-J 1 



^->3e, 
^ Ti-»e Ti, 
p Au->e Au excluded 
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Al->e Al, 
/i Ti->e Ti proposed 
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<1 
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^ Al->e Al, 
/; Ti-»e Ti proposed 



(f) 



FIG. 3: (a)-(d): Parameter space excluded (projected) by various bounds (expecting sensitivities) 
on \i LFV processes through contributions from photonic penguins. Note that solid, dashed, dot- 
dashed and short-dashed lines denote results from /j, — > e"f, /i — > 3e, /uTi — > eTi and /iAu(Al) — > 
eAu(Al) processes, respectively, (e) and (f): Same as (a)-(d), but through contributions from 
Z-penguins. 



(ii) For m^/m^ > 0.1, the allowed bands raise with the mass ratio. The muon g — 2 is more sensitive 
to the diagram with ip interacting with a photon [as depited in Fig. [T|(h>)] than to the other diagram, 
hence, the constraint on —Q^Re{g*^ R g^i)/m^, is severer than the one on +Q ( pRe(gJ lR g^L)/m^. 
Indeed in the large m^/m^ region, the suppressions from a large mass should be larger in 
diagrams with <p interacting with a photon and, hence, they require larger couplings to compensate 

11 



/j->3e, 
ji Ti->e Ti, 
ju Au-»e Au excluded 



fi->3e, fi Al->e Al, 
ju Ti->e Ti proposed 



500 

m^GeV) 
(a) 



/i->3e, 
Ti->e Ti, 
p Au->e Au excluded 



^->3e, fi Al->e Al, 
\i Ti->e Ti proposed 



100 150 200 300 500 700 1000 



m (GeV) 
(c) 



/j->3e, 
fi Ti->e Ti, 
fi Au-»e Au excluded 



/i->3e, ii Al-»e Al, 
fi Ti->e Ti proposed 



m (GeV) 
(b) 



/j->3e, 
fi Ti-»e Ti, 
fi Au->e Au excluded 



ju->3e, fi Al->e Al, 
fi Ti->e Ti proposed 



m/GeV) 
(d) 



FIG. 4: (a)-(d): Parameter space excluded (projected) by various bounds (expecting sensitivities) 
on \i LFV processes through contributions from Z penguins with different choices of mass ratios. 
Note that dashed, dot-dashed and short-dashed lines denote results from fi — > 3e, //Ti — > eTi and 
/iAu(Al) — > eAu(Al) processes, respectively. Note that these plots also apply to the R «-)■ L cases. 



the effect. We see that the mass ratio m^/m^ can go up to 200 (1000) along the 4> (VO band, which 
corresponds to allowing to be as large as 20 (100) TeV in the extreme situation. 

As noted previously, in Fig. [2^b) we also show the excluded region of iQ^^Im^g^g^l/m^ 
from the muon EDM bound. We see that the bound is three order of magnitude higher that the 
allowed ±Q (f,^Re(g^ R g^L) / bands. To constrain the former to the level of the latter, the EDM 
sensitivity needs to be improved. In fact, some proposed EDM searches (see, for example |llj ) are 
aim at a 4 to 5 order of magnitude improvement on the sensitivity and may be able to probe the 
imaginary part of g^gpL better than its real part. 

We now turn to /i LFV processes, including /x — > e'y, /x — > 3e, /iTi — > eTi and /iAu(Al) — > 
eAu(Al) transitions. In Fig. [3] and [4j we show the parameter space excluded by various bounds 
and the one corresponding to projections from the expected sensitivities on these /t LFV processes, 
through contributions from photonic and Z penguin diagrams. To be specific, for the proposed 
sensitivities, the conservative values of the future sensitivities quoted in Table 1 are used. Note 
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(c) (d) 
FIG. 5: Ratios of rates contributed from photonic penguins. 



that the photonic penguins contribute to [i — > 3e and /xiV — > ciV through the so-called F% penguin, 
which is similar to those contributing to Aa M and fj, — > ej, and the Fi penguin, while the Z- 
penguins only contribute to fj, — >• 3e and /uiV — > eN decays. Note that the Z-penguin amplitudes 
contribute through the \g,j,R(L)9eR(L)AT 3i> \ and \g^R(L)9eR(L)^R(L)\ parts. The former contribution 
is a function of the mass ratio m^/rn^, while the latter one depends on both (fi and ip masses. The 
resulting constraints are plotted in Fig. [3^e), (f) and[4j It should be noted that Fig. [3] and [4] can 
still be useful if the experimental bounds change. For example, if a bound is reduced by a factor 
of k, the new plot can be easily updated by reducing the present plot by a factor of yk. 

Note that the combinations of couplings \Q<t,^g^R{L)9eR(L)\, \Q<j>,^9^R(L)9eL(R)\ (from photonic 
penguin) and \9^R(L)9eR{L)^T^\, \g^R{L)9eR{L)i^L{R)\ (from Z penguin) have different sensitivities 
on the experimental constraints and the sensitivities change with m^^. For heavier Jfuy,, the 
contributions from the photonic penguins decrease and the Z-penguin contributions, where the 
non-decoupling effect is working, dominate. By comparing Fig. [3^a), (b) to Fig. [3](e), (f) and[4j 4 
we find that (i) in the range of ^ 0(100) GeV the photonic penguin contributions dominate over 



4 Note that the plotted quantities in these figures have different powers of . 
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FIG. 6: Constraints on parameters which contribute through box diagrams to the ^ + — > 3e process. 
Solid lines denote the fj, — > 3e constraint or expectation in the Dirac case, while the dashed lines 
denote the Majorana case. Note that these plots also apply to the R L cases. 



the Z-penguin ones, (ii) in the range of the 0(100) GeV < < 0(100) TeV, the photonic penguin 
contributions from the Q(f>^g^R{L)9eL(R) terms dominates over the Z penguin contributions, which 
are, however, still larger than the photonic penguin contributions from the Q(f>^g^R{L)geR(L) part, 
and for (iii) 0(100) TeV the Z penguin contributions dominate. The role which Z penguin 

plays is emphasized in [12]. 

Since NP contributions to Aa^ and the fi + — > e + 7 decay are from similar diagrams, it will be 
useful to compare them. Using Fig.[2^b),[3](c) and[3^d), the present data on Aa M and B{ji + — > e + rf 
lead to 



9fj,R(L)9eL(R) _ 9eL(R) 



< 6.1 x 10" 5 ~ A e 



(35) 



9fiR9fiL 9fj,L(R) 

where we define A = 0.2. This ratio is much smaller than any known coupling ratio and mixing 
angle among the first and second generations. For example, the mass ratio m e /m^ ~ A 3 ~ 4 , quark 
mixing in CKM matirx V u d = sm.6 c ~ A, neutrino mixing sin #^12 ~ y/X are all larger than the 
estimated g e L(R)l 9ij,l(r) coupling ratio. It seems that the present case is unnatural. 

We see from Fig. [3^a)-(d) that the present bound from fi — > &y surpasses all other bounds. In 
particular, even the parameter space to be probed by the proposed [i — > 3e sensitivity is mostly 
excluded by the present fj, — > erj bound. This can be understood by using Fig. [5j where rate 
ratios of various modes through photonic penguins are given. We see that the ratios of LFV rates 
with respect to the (i —> cy rate are all less than unity. Furthermore, we recall that the present 
experimental bounds on LFV rates are of similar orders of magnitudes (see Table 1). Therefore, 
the present bound on the /i — > ej rate provides the most severe constrain. 

Taking a closer look at Fig. [5j we see that, from Fig. [5] (a) and (b), the \Qip,^g^R(L)9eR(L)\ terms 
give B(fi ->• erf > B((J,N eN) > B(fi -)• 3e) and, from Fig. [5] (c) and (d), the \Q4>^g^R{L)g e L(R)\ 
terms give B(fi 3e)/B(/z erf ~ 0.006, 0(//Ti eTi) /B((i erf ~ 0.004 and B(/xAl 
eAl) /B(fj> — > erf ~ 0.003, where the first ratio is consistent with Ref. [6]. Hence, for g^R{L)geL(R) 
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dominating models, the latest MEG bound implies 

B(li+ -»• e+e+e") ~ 0.006 x B(p + -»• e+7) < 1.4 x 10" 14 , 

B(fiN -> e2V) ~ O(10- 3 ) x 2%+ -> e+ 7 ) < O(10" 15 ), (36) 

for Af =Au, Al and Ti. The above expecting limits are about two to three orders of magnitudes 
below the present experimental sensitivities (see Table I) and make the searches on LFV in the 
muon sector challenging in this case. 

As noted in Sec. II A, it is possible to have box diagrams with the electron line in Fig. [TJc), (d) 
replaced by a quark one and contribute to fj,N — > eN conversions in some cases. The correlation to 
the (jl — > 3e rate will be modified. We will discuss more on this situation in the discussion section. 

Note that Z-penguins give different rate ratios (not shown in Fig. [5]), with £>(/xAl — > eAl)/£>(/x — > 
3e) ~ 10, B(/iTi -> eTi)/B(n -> 3e) ~ 20 and £(/xAu -> eAu)/B(/i -> 3e) ~ 40 roughly indepen- 
dent of the masses m^. This pattern is different from the photonic penguin case as shown Fig.[5j 
These rate ratios will be useful for identifying the underlying NP contributions. 

In Fig. [6j we show the constraints on parameters which contribute through box diagrams, as 
depicted in Fig. [I] (c) and (d), to the fi + — > 3e process. Both Dirac and Majorana cases are shown. 
We see in Fig. [6] (a) that there is cancellation in the Majorana case and the sensitivity on the 
parameters is relaxed. 

Note that constraints on the same combinations of parameters can be obtained from penguin 
processes, including fi — > erf, Aa e , EDM, and the purturbative bounds, as well. They are also 
shown in Fig. [6j We see that these constraints are usually much stronger than the ones from the 
box diagrams, except for \g^,R(L)9 e R(L)9eR(L)9eR(L)\/i , n^ in the low m^/m^ region. In particular, 

the fi -> e7, Aa e and the electron EDM constrain (Q 9 i,v>|9 At /?(L)S' e L(i?)l/wv)( ( 3</»,'/'fi , e^(L)5eL(i?)l/"iv) 
much deeper than \g flR { L )9eR(L)9eL(R)9eL(R)\/mf l} from the box diagrams. Hence, in general, these 
box diagrams do not play a major role in the fi + — > 3e decay. 



B. Case II 

We now turn to case II. In Fig. [7J the allowed regions for = fQ(p,^9tiR9tiL^(^RL)/i^/'>TT'ip con- 
strained by the measured Aa M , with exclusions of \9^l9^rSrl\ > 4tt and < 100 GeV, are 
shown. Excluded and projected parameter space of \Q ( ) > ,' l j}g l j,Rg l j,L^ 1 {^RL)n^\/'mf from the muon 
EDM bound and the expected sensitivity are also given on the same plot. For the plots of the 
allowed regions for ±Q(j,^\g^L(R)\ 2 l^tfpi one is referred to Fig. |2|a), as they are common in both 
cases. . 

Comparing Fig. [7] to Fig. [2](b), we see that the allowed parameters in the upper m^/m^ region 
are similar. In contrast, they are relaxed substantially in the lower m^/m^ region in the present 
case. To reproduce the measured Aa M , we need to have 0(1O~ 2 ) < m^/m^ < O(10 2,3 ), where the 
minimum of the mass ratio is much higher than the one in the previous case. Recall that in case 
I, for m^/m^ < 0.1, the allowed parameter region for ±Q ( p^Re(g^ R g tl L)/mjp are horizontal bands 
around 4 x 10~ 6 GeV -1 and m^/m^ can be as low as O(10" 5 ). From Fig.g we see that as we move 
downward along the m^/m^p axis, the bands for the allowed regions for TQ^^g^Rg^L^iSR^^/m^ 
bend upward in the low mass ratio region (m^/m^ < 1) and the above parameters can be as large 
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FIG. 7: Allowed parameter space for TQ^^g^Rg^L^^R^fifi/^nip constrained by Aa M with ex- 
clusion of \g^L9fiR^Ftx\ > 47r and m^^ > 100 GeV. Excluded parameter space (shaded regions with 
solid or dashed lines) of |Q^ i ^<7 /J _R5' Ai Llm(5^) A j /J /m^| from the muon EDM bound is also shown. 



as 10 -3 GeV^ 1 , which is three orders of magnitude higher than those in case I. We also note that 
the mass ratio m^/m^ cannot be smaller than 10 -2 as the bands quickly run into the shaded 
rigions, which correspond to the excluded < 100 GeV and \g ^Rg At iRe((5 >rl) w | > Air regions. 
In the present case, the mass of ip cannot be larger than a few tens TeV, while in case I it can 
be as high as few thousand TeV in the extreme situation. The built-in cancellation mechanism 
reduces the amplitudes effectively and a too heavy ip is incapable to produce a large enough Aa M . 
The effect of the cancellation is important in the low m^/m^ region and, consequently, relaxes the 
constraints on parameters. In fact, we expect to see the very feature in other penguin contributing 
channels as well. 

In Fig.[8j we show the constrained and projected parameter space through penguin contributions 
by considering the experimental bounds and the proposed sensitivities. 5 We note that the photonic 
penguin contributions via the 5rl term dominate over other contributions for below 0(1O 3 ) 
TeV. For mw, beyond that the Z-penguin contribution takes over. However, from the previous 
discussion on the muon anomalous magnetic moment, we see that to account for the measured Aa^, 
m^p cannot be heavier than few tens TeV. Hence, the Z-penguin contribution will be subdominant 
in this case. 

By comparing the constraints from Aa M and B{(jl —> cy) as shown in Fig. [7J [8]^c) and|8^d), we 
obtain 

g^R{L)9eL(R)^[{^RL{LR))^e] _ g e L(R) ^[(^RL(LR))fie\ 



g^Rg^L^[(8 RL )^} 



< 4.2 x 10~° ~ A 



6 



(37) 



If we estimate g e UR)l 'g^LfR) by using the lepton mass ratio m e /m^ 



X 



3~4 



we see that a mixing 



angle ratio of Re[(<5^(i^)) Ate ]/Re[(5ijL) AtAt ] ^$ A 2 ~ 3 , which is not unnatural, can easily satisfy the 
above bound. In this respect, case II is more reasonable and natural than case I, where the coupling 



ratio is highly hierarchical [see, Eq. (35)]. 



For Z penguin contributions, only those from AT^^, are shown, since the Hl(r) ones are highly suppressed. 
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FIG. 8: Same as Fig. [3j but now in case II. 



It is interesting to see from Fig. [8] that the fj, — >■ erf bound is not always the most stringent one. 
The bounds from /i — > 3e and /iiV — > eN in Fig. [8^a) are almost the same as those in case I [see 
Fig. |3^a)], but the bound from [i — > ej is relaxed up to more than two orders of magnitudes in 
the low m^/m^ region and becomes less severe than other bounds. Similarly, comparing Fig. [8|c) 
with Fig. [3]^c), we see that in the low m^/m^ region both bounds from [i — > e^ and fi — > 3e are 
relaxed up to three orders of magnitudes, while the changes on those from fiN — > eN are mild. We 
can infer that, similar to the Aa^ case, the F2 (photonic) penguin amplitudes exhibit cancellations 
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FIG. 9: Same as Fig. [5} but in case II. 



in amplitudes in the low m^/m^ region and relax the constraints from [i — > ey significantly, while 
the cancellations in the Fx penguin contributions in fi — > 3e and (jlN — > eN processes are mild. As 
a result the bounds from liN — > eN approach the — > ey bound in this case, while in the previous 
case these two bounds are always apart. 

The ratios of photonic penguin contributing rates plotted in Fig. [9j show that B(fj,N — > 
eN)/B(fi — > ey) and BQu — > 3e)/B(fi — > ey) are enhanced compared with those in Fig. [5} 
In Fig. |9^a) we see that the ratios can be enhanced up to three orders of magnitudes, in 
Fig. [9^c) the B(fiN — > eN)/B(^L — > ey) ratio is enhanced by one order of magnitude, while the 
B(n — > 3e)/B(fi — > ey) ratio does not change much. It is very interesting that the rate ratio 
B(fxN — > eN)/B(fi — > ey) from the ff^R(L)5eL(R) term is enhanced and different from case I. 

We see in Fig. [8] that parameters with 5rl [as shown in (c) and (d)] are most constrained by 
data. It is likely that these parameters give dominate contributions to LFV processes. Using 
Fig. [9^c) and (d) we find that the present bound on — > ey allows 

B(jtN -»• eN) < 1(T 13 , (38) 

which is close to the present bounds (see Table [T]). Therefore, the search on these processes could 
be very interesting. 
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In Fig. 10, we show the constraints on parameters which contribute through box diagrams to 
the /i + —> 3e process in this case. Although we also see some relaxations on parameters, the main 
conclusion remains similar to that in case I. 



IV. DISCUSSIONS 

A. Flavor violating Z-decays 

Lepton flavor violating Z — > fi^e^ 1 decays are highly related to \i — > 3e, [i — > and fj,N — > eN 
processes via the Z penguin contributions. The Z — > I' I decay rate is given by 



rrir. 



Z \2\ 



nz^l'l) = ^{\gf jlL \ l + \9R'R 



(39) 



where the dimensionful coefficient g%, M is the same one used in Eq. ( ]12j ). Using Eq. (p9J) and 
the results in the previous section, we find that the present bounds from — > 3e, /iTi — > eTi and 
^Au — > eAu processes constrain 



B(Z -> /x T e ± ) < 4 x 10~ ia , 7 x l(T i4 , 6 x 10 



13 



-14 



-15 



(40) 



respectively. Note that the above equation holds in both case I and II. In any case, these constraints 
are far below the present limit, B(Z — > fJ^e^ 1 ) < 1.7 x 10~ 6 . 



B. Box diagrams involving quarks 

Form the explicit assignment of gauge quantum numbers of tp and 4> as shown in Appendix lAl 



we see that it is possible to have 4> couples to quarks [see Eq. (A5)]. These interaction can generate 
additional contributions to fj,N — > eN conversion precesses through box diagrams similar to those 
in the fi — >• 3e ones as shown in Fig. 1(c) and (d), but with the (lower) electron line replaced by a 
quark line. 
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To have interaction with quarks, only rather specific choices of ip and gauge quantum numbers 
are allowed (see Appendix [A|) . In particular, the case of i^r : (1,1,1) and (j)L '■ (1,2,-1/2) are of 
interest, for the fermion field is a SM singlet. In below we will use this case to illustrate the 
contributions from the additional box diagrams. 

The interacting lagrangian in this case is 

£int = giL^RLLi4>Li + 9uQLiUR<t>L t + 9dQ "hi^RUj^hj + h.C, (41) 

where Ql and Ll are the quark and lepton doublets, respectively. Note that only the lower 
components of Ll and (pi are relevant to this analysis. The box diagrams give 

9R V (q) = 0, g L v(d) = 0, (42) 

and 

9Lv(u) = - \ G ( m l'°' m % m l)(9*eL9^L)(9t9u +0d0d)}> ( 43 ) 

where quark masses have been neglected. Note that the box diagrams also give the so-called 
9lp{u) term, which, however, does not contribute to conversion rates [13]. The resulting fiN — > eN 



conversion rates can be calculated using Eq. (18). 



Experimental limits on conversion rates are used to constrain couplings and masses. The result 



is shown in Fig. 11 The correlation between [iN — > eN conversions and the /U + — > 3e decay are 
lost. In fact, we see that the constraints on \gfj,L9eL\{\9u\ 2 + \dd\ 2 ) /m^ from present limits on \x — > e 
conversion rates are similar to the constraints on \g l _ L L9eL\\9eL9eL\/ m2 b from the /x + — > 3e bound 
(see Fig. [6]). Therefore, we may be able to see fiN — > eN conversions sooner than the fi + — > 3e 
decay, if g u d is larger than g e i, and vice versa. The fiN — > eN conversion rates need not be highly 



suppressed as noted in Sec. Ill A 
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C. Some other cases 

A similar analysis can be preformed by replacing the spin-0 particle by a spin-1 one in the loops. 
It will be interesting to compare it to the present work. However, gauge invariant and triplet vector 
couplings will complicate the analysis. The study will be given else where. 

We expect to find results similar to case II, but with cancellation at work in low region, if 
we introduce the built-in cancellation mechanism in the ip sector, instead of in the <ft sector. 



V. CONCLUSIONS 

In conclusion, we use a model independent approach in this analysis, where these processes are 
considered to be loop-induced by exchanging spin-1/2 and spin-0 particles. We explore two com- 
plementary cases, which has no or has an internal (built-in) cancellation mechanism in amplitudes. 
Our main results are as follows: 

(a) Bounds from rates are used to constrain parameters, such as coupling constants and masses. 
These constraints can be easily updated by simple scalings, if the experimental situations 
change. 

(b) The muon g — 2 data favors non-chiral interactions. 

(c) In |U + — > e + e~e~ and [i~ N — > e~N processes, the Z-penguin diagrams may play some role, 
while the box diagrams contributions to the fi + — > e + e~e~ rate are usual highly constrained. 

(d) Z-penguin contributions can be constrained from — > e + 7 and fj,~N — > e~ N bounds. It 
can then be used to constrain the Z — > e^fj, rate by 7 to 8 orders of magnitudes lower than 
the present experimental bound. 

(d) In the first case (without any built-in cancellation mechanism), using the recent fi + —> e + 7 
bound, we find that [i + —> e + e~e~ and fi~N — > e~N rates are bounded below the present 
experimental limits by two to three orders of magnitudes in general. In some cases, the 
above expectation on low fi~N — > e~N rates can be relaxed, as additional box diagrams 
involving quarks contribute to fj,~N — > e~N processes. 

(e) Furthermore, by comparing Aa M and B(/j, — > ej) data, the couplings of and g e are found 



to be highly hierarchical [see Eq. (35)]. Additional suppression mechanism should be called 
for. 

(f) In the second case (with a built-in cancellation mechanism), mixing angles can provide 
additional suppression factors to satisfy the Aa M and B(fi — > cy) bounds without relay only 
on highly hierarchical g e and g^ couplings. 

(g) In addition, although the /x + — > e + e~e~ rate remains suppressed, the bounds on — > 
e~N rates, implicated from the MEG /i + —> e +r y bound, can be relaxed significantly in the 
second case and can be just below the present experimental limits. 



21 



TABLE II: Some examples for the assignment of the quantum numbers of ipL,R and 4>l,r- 





4>L 


4>L 






(1,1,^) 


(l,2,Y R -l) 


(1,1, ^l) 


(i,i,y L - 


i) 


(1,2,**) 


(1,1,1^-1) 


(i,2,y L ) 


(i,2,y L - 


i) 


(3(3),i,y fl ) 


(3(3),2,y B -|) 


(3(3),i,y L ) 


(3(3),i,y L 


-i) 


(3(3),2,y fl ) 


(3(3),i,y fl -§) 


(3(3), 2, y L ) 


(3(3), 2, y L 


-i) 
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Appendix A: Gauge quantum numbers of (p and ip 

The tp — 4> — I lagrangian, 

Ant = g' L (^R4>* L )i(L L )i + g' R tpL^ R lR + h.c, (Al) 

where i is the weak isospin index, is gauge invariant under the SM gauge transformation. Recall 
that the lepton quantum numbers under SU(3) xSU(2) xU(l) are given by 

Lz:(l,2,-i), Ir: (1,1,-1). (A2) 

The gauge invariant requirement implies that we must have the following quantum number assign- 
ments for these combinations: 

«L:(l,2,i), ^ : (l,l,l). (A3) 
Consequently, the gauge quantum numbers of ip and (j) are related as following: 

i/> R : (c R , 2I R + 1, Y R ), <p L ■ (cr, 2{I R ± -) + 1,Y R - -), 

tfj L : (c L ,2I L + l,Y L ), </> r :(c l ,2I l + 1,Yr-1). (A4) 

Some examples for the assignments of the quantum numbers of ipL,R an d <Pl,r are given in Table |Tlj 
Note that in the cases of Ir = 0, Y R = and II = 1/2, Yl = 1/2, <pL and 4>r can couple to 
quarks, respectively, through 

QuUR&Li, QudR^m, QLiURCij^Rj, QL t dReij(jJ Lj , (A5) 

where ey is the antisymmetric tensor. It is easy to see that the above terms are indeed gauge 
invariant by using QlQr ■ (1,2, —1/6 + Q q ) and Eq. (A4). 



Appendix B: Loop functions and input parameters 

The loop functions used in this work are defined as 

1 

12(o - b) 



Fi(o,6) = 1 ( 2a 3 + 3a 2 b - Qab 2 + b 3 + 6a 2 b In - 
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TABLE III: Parameters of overlap integrates and total capture rates w cap t taken from |X3|, [T4"] . 



D(m 5 J 2 ) 



5/2n 



m 



5/2^ 



Wcapt(10 6 S" 



!3 7 A1 

iTi 

205 Tl 
81 11 



0.0362 
0.0864 
0.189 
0.161 



0.0161 
0.0396 
0.0974 
0.0834 



0.0173 
0.0468 
0.146 
0.128 



0.7054 
2.59 
13.07 
13.90 



F 2 (a,b) 
F 3 (a,b) 
Gt{a,b) 
G 2 (a,b) 
G 3 (a,b) 
F z (a 1 ,a 2 ,b,b, c) 

F z (a,a,bi,b 2 ,c) 

G z (ai,a 2 ,b) 
F(a, b, c, d) 



G(a, b, c, d) 



2(a- 


6) 3 \ 


1 




2(a- 


6) 3 ( 


1 




36(a- 




1 




36(a- 




1 





3a 2 + 4a6 - b 2 - 2a 2 In - 

a 

b 



a 2 -6 2 + 2a61n 



a 



36(o - 6) 5 



(a - 6)(lla 2 - lab + 26 2 ) - 6a 3 In 



(a - b) (16a 2 - 29ab + 7b 2 ) - 6a 2 (2a - 36) In 

(a - 6)(17a 2 + 8a6 - 6 2 ) - 6a 2 (a + 36) In ^ 

a) 



a\(2Ja\a 2 - a x ) a x a 2 (2Jaia 2 - a 2 ) a 2 
In h — - ~ — m 



2(«i — a 2 )(ai — b) c 2(a\ - a 2 )(a 2 - 6) c 

b(2Ja\a 2 - 6) 6 
In 



2(ai — 6)(a 2 — 6) c 

,2 



3 ^ In" ^ ln^ 1 

~4 + 2(a-6i)(a-6 2 ) n c ~ 2(a-6i)(6i - 6 2 ) " 7 

, b 2 62 

2(a-6 2 )(6i -6 2 ) n c ' 



aiyja\a 2 



In 



f/i 



«2 V /fl l a 2 



hi 



«2 



(ai — a 2 )(ai — b) b («i - a 2 )(o 2 - 6) 6' 



6ya6 6 cvo6 c 

In In — 

(a — 6) (6 — c)(6 — d) a (a — c)(b — c)(c — d) a 



d\ ab , d 

-I In - 

(a — d)(b — d)(c — d) a' 

6 2 6 c 2 , c 

In — I In — 

(a — 6) (6 — c)(6 — d) a (a — c)(b — c)(c — d) a 

d 2 , d 

m 



fBl) 



(a — d)(b — d)(c — d) a 

Note that these loop functions are dimensionful and the dimension of G3 is different from others. 
We do not have the expression of Fzifli, a 2 , 61, 6 2 , c), since in Sec. II C only a\ = a 2 = a and/or 
61 = 6 2 = 6 are needed. Both expressions of F z give identical result in the a\ = a 2 = a and 
61 = 6 2 = 6 case. 



The numerical values of D, V and u; capt used in Eq. (18) are collected in Table 
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